Algebraic spin liquid in an exactly solvable spin model 
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We have proposed an exactly solvable quantum spin-3/2 model on a square lattice. Its ground 
state is a quantum spin liquid with a half integer spin per unit cell. The fermionic excitations are 
gapless with a linear dispersion, while the topological "vison" excitations are gapped. Moreover, 
the massless Dirac fermions are stable. Thus, this model is, to the best of our knowledge, the first 
exactly solvable model of half-integer spins whose ground state is an "algebraic spin liquid." 



The term "spin liquid" is widely used to give sharp 
meaning to the more general intuitive notion of a Mott 
insulator; a spin liquid is an insulating state that cannot 
be adiabatically connected to a band insulator, i.e. to 
an insulating Slater determinant state. In a system that 
preserves time reversal symmetry, any insulating state 
with an odd number of electrons (or a half-integer spin) 
per unit cell is a spin liquid. Interesting proposals [l|, l2| 
have been made concerning the relevance of such states to 
the theory of high temperature superconductivity in the 
cupratcs and other materials. Indeed, various spin-liquid 
phases have been proposed, which arc distinguished by 
the character of any gapless spinons and the exchange 
statistics of the topological "vison" excitations. 

Since they are new and "exotic" quantum phases of 
matter, it is desirable to construct solvable models with 
short range interactions with stable spin-liquid ground- 
state phases. A breakthrough occurred when Mocssncr 
and Sondhi Q demonstrated the existence of a gapped 
spin liquid ground state in the quantum dimer model [H , 
analogous to the short range version of the RVB state 
0,0]. An exactly solvable spin-1/2 model in a gapped Z 2 
spin liquid phase was later constructed by Wen [7|] . How- 
ever, much of the recent interest, spurred in part by the 
possible observation of such a state in re-(ET) 2 Cu2(CN)3 
H, 0, G3 and Zn(Cu) 3 (OH) 6 Cl 2 HH, has focussed on 
spin- liquids with gapless spin excitations, so-called "al- 
gebraic spin liquids." 

The exactly solvable Kitaev model on the honeycomb 
lattice [HI can exhibit gapless excitations. However, be- 
cause the honeycomb lattice has two sites per unit cell, 
this model has an integer spin, hence an even number of 
electrons per unit cell. In the present paper we construct 
an exactly solvable model, in much the same spirit as the 
Kitaev model, whose ground state is a spin liquid with 
an odd number electrons per unit cell and stable gapless 
fermionic excitations. To the best of our knowledge, this 
is the first exactly solvable model with this sort of spin 
liquid ground state -algebraic spin liquid. 

The Kitaev model has a spin-1/2 on each site of a 
trivalent lattice, where the coordination number is dic- 
tated by the existence of three Pauli matrices. In order to 
study a model on a square lattice, we instead consider a 
model with a spin-3/2 on each lattice site. The resulting 



larger Hilbert space, with 4 spin polarizations per site, 
permits us to express the model in terms of the 4x4 an- 
ticommuting Gamma matrices, T a (a — 1, • ■ • , 5) which 
form Clifford algebra, {r a ,r 6 } = 26 ab . Specifically, the 
5 Gamma matrices can be represented [l3} by symmetric 
bilinear combinations of the components of a spin 3/2 
operator, S a , as 
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Model Hamiltonian: We define our model on a square 
lattice, with a spin-3/2 on each site, and corresponding 
r matrices defined as in Eq. fTJ). In terms of these, 
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where r afc = [r a , T b ]/(2i) and i labels the lattice site 
at Yi = (xi,yi). We call this model the Gamma ma- 
trix model (GMM). Suppose that the square lattice has 
N = L x L y sites, where L x and L y are the lattice's linear 
sizes and are assumed, for simplicity, to be even in this 
paper. Moreover, we consider periodic boundary condi- 
tions. Obviously, the GMM can be written explicitly as a 
spin-3/2 model. The GMM model respects translational 
symmetry and time reversal symmetry (TRS). It does not 
have global spin SU(2) or even U(l) rotational symmetry, 
but is invariant under 180° rotations about the z axis in 
spin space, i.e. it has global Ising symmetry. The lack of 
any continuous symmetry implies that the fermionic ex- 
citations (discussed below) do not have well-defined spin 
quantum number, so their relation to "spinons" is un- 
clear. A feature of the model which makes it solvable is 
an infinite set of conserved "fluxes:" [W,-,?i] = for any 
i and [W^Wj] = 0, where W, = r 13 ^ 23 
the plaquette operator at site i. 

Note that, in contrast to the spin-1/2 Kitaev model on 
the honeycomb lattice, the GMM has an odd number (3) 
of electrons per unit cell. The present model in the limit 
J x = J y = J' x = J' y and J5 = is similar to a model pro- 
posed by Wen in Ref. 14j, where, however, the Gamma 
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matrices were constructed from two spin- 1/2 operators 
on each site, and therefore behave differently under time 
reversal than in the present realization. Moreover, the 
present model generically does not possess the special 
symmetries of Wen's model that are responsible for some 
of the behaviors of that model. 

Fermionic Representation: Spin-3/2 operators can be 
expressed as bilinear forms involving three flavors of 
fermion operators, S z = a^a + 2b^b — 3/2 and S + = 
VSf^a + V3af + 2aJb, subject to the constraint that the 
physical states are only those with odd "fermion parity," 
(-1)* = -1, where N = f^f+ata+tfb. Rather than us- 
ing this representation in terms of three Dirac fermions, 
we will directly represent the Gamma matrices in term 
of a related set of 6 Majorana fermions: 
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where , di, and d\ are Majorana fermions on site i. 
Six Majorana fermions form an eight dimensional Hilbcrt 
space, which is an enlarged one from the physical Hilbcrt 
space of an spin-3/2. In terms of spin-3/2 operators, 
rirfrfrfrf = -1 for all i. Consequently, all allowed 
physical states in term of Majorana fermions must 
satisfy the following constraint, for all z, 
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In terms of Majorana fermions, the Hamiltonian in the 
enlarged Hilbert space can be written as 



H = 



+J' x u i>x id' i d' i+X + JyU^yid'id 



Jbidid 1 ; 



(5) 



where u^ x = —ic}c 2 + ~ and M^y = —icfc* + ~. It is obvious 
that Ui t \ are conserved quantities with eigenvalues m^a = 
±1, A = x, y. The enlarged Hilbert space can be divided 
into sectors {u}. In each sector, the Hamiltonian Eq. ([5|) 
describes free Majorana fermions: 
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where Ui t \ are emergent Z 2 gauge fields. The Z 2 gauge 
transformations arc given by di — > A^d,-, d[ 



Aid' i: and 



AiUi,\Ai + \, where Aj = ±1. In the enlarged 
Hilbert space, the eigenstates = \ip} c ® \ip)d,d' can 
be written as a direct product of a part that involves 
the fermions and a part that involves the d and d' 
fermions, respectively. Here \ip} c is defined by w^aIV^c = 
Ui,\\ip) c and \ip)d,d' are eigenstates of Eq. ([6|) with the 
corresponding u.i t \'s. 

The spectrum of Ti({u\) depends only on gauge invari- 
ant quantities - the flux on local plaqucttcs, cxp(i</>j) = 



and two global fluxes exp(i(j> x ) = 
Hi( yi =i) and exp(i(j> y ) = rii(a ii =i) Wi'.j/j wncr G and 
<l>x,y = 0, 7r. [Note that the previously defined Wi = 
— exp(i(j>i) .] It is obvious that J2i 4>i = (mod 2tt), so 
there are N — 1 independent local fluxes. Including the 
two global fluxes, the number of independent fluxes is 
N + 1. Since there are 2N Z 2 gauge fields, the num- 
ber of different gauge field choices corresponding to each 
flux sector {(/)} is 2 N ~ 1 . In other words, in the enlarged 
Hilbcrt space, each state is 2 Af_1 -fold degenerate. Note 
that in the thermodynamic limit, the energy is indepen- 
dent of the two global fluxes, which gives rise to a fourfold 
topological degeneracy of the physical ground states. 

Projection operators: Most of the states in the en- 
larged Hilbert space are not physical states. To obtain 
a physical eigenstate, we must find a linear combination 
of the degenerate eigenstates which is simultaneously an 
eigenstate of every Di with eigenvalue 1 . This is realized 
through the projection operator P: 



i*}=^}^n[( i+ ^)/ 2 
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Clearly, Eq. ^ implies A|^) = |^) for any i. Explicitly, 
P is given by 
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Di acting on a direct product state |?/>) is equivalent 
to a gauge transformation on site i. A subtlety here 
is that there are 2 N operators in the sum in Eq. l[8]). 
but only 2 Af_1 incquivalcnt gauge transformations. In 
fact, D = Y[ i Di implements a gauge transformation on 
every site, thus leaving all gauge fields unchanged. It fol- 
lows that P = P'{\ + D), where P' includes all inequiv- 
alent transformations. Since [D,H] = and D 2 = 1, 
D\ip) = ±\ip). Moreover, D = Hi IL 
In terms of the number Nf of the Dirac fermions, 



fj^i^dj+id'A/2, N f 



(9) 



and fluxes, N^,, we obtain D = (—l) N * , (—l) Nf . Nf is con- 
served modulo 2 by the Hamiltonian. N^ is defined by 
dividing the plaqucttes into two sublattices and counting 
the number of 7r-fluxes through one or the other sublat- 
tice Q. 

Thus, depending on the fermion and flux parity, P ei- 
ther annihilates a given direct product state, \ip) , or 
maps it to the equal weight linear superposition of all 
gauge transformations acting on \tp). For instance, when 
there is a 7r-flux through each plaquette, which, as dis- 
cussed below, is the ground state sector, D = (—l) N f, 
i.e. all physical states must have an even number of 
fermions. Conversely, where 1 7r-flux is added to each 
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sublattice, only states with Nf = odd survive projec- 
tion. Conserving the parity of fermion number reflects 
the fact that physical fermionic excitations are created 
by non-local (string) operators [HI, . 

TT-flux state and gapped visons: In each flux sector 
{fa}, the lowest energy of the Hamiltonian is denoted 
by E ({fa}). The ground state energy of the model is 
achieved by minimizing Eo({fa}) with respect to {fa}. 
Formally, by integrating out the fermions, an effective 
action for the Z 2 gauge fields can be derived. However, 
in general, it is nontrivial to obtain an explicit form of 
the effective action. 

For J 5 = 0, fortunately, there is a theorem due to Lieb 
[IH which implies that the energy minimizing flux sector 
of a half-filled band of electrons hopping on a planar, 
bipartite lattice is tt per square plaquctte. We define 
this uniform 7r-flux, the ground state sector, as being 
vortex-free. Z 2 vortex excitations, "visons" , are defined 
to be plaquettes with fa = 0. Due to the constraint 
Y] ■ fa = (mod 2tt), only an even number of visons 
are allowed. It is special for this model that visons do 
not have dynamics even though they interact with one 
another. (Numerical calculations of the two vison energy 
reveals that the interaction between visons is repulsive.) 
The minimum energy cost of creating two visons in the 
ground state is always finite and is defined as the vison 
gap energy A.„ ~ (^/| J x J y \ + ^\J' x J' y \). Since the visons 

are gapped, the uniform 7r-flux is still the ground state 
sector as long as J5 is much smaller than A„ . In the rest 
of paper, we restrict our attention to cases (e.g. small 
J5) in which the ground state lie in the 7r-flux sector. 

Spin correlation: The state with uniform 7r-fluxes pre- 
serves the translational symmetry and TRS of the model. 
To prove the ground state of the model is a spin liquid, 
we need to show that there is no magnetic order and 
that the spin correlations are short-ranged. Spin-3/2 op- 
erators can be expressed as bilinear forms of Majorana 
fermions. For instance, Sf = ic\c\ + \ic\c 2 . The action 
of Sf on the ground state Y$>) creates visons in the four 
surrounding plaquettes around site i [l(|- The effect of 
Sf and Sf on the ground state is similar. Because visons 
are non-dynamical in the present model, (^!\SfSj\^f) is 
identically zero unless sites i and j are nearest neighbors, 
i.e. the spin correlations are unphysically short-ranged. 
Presumably, if additional small, local terms are added to 
the Hamiltonian, perturbative corrections to this corre- 
lator would lead to a finite correlation length. 

Gapless fermions: To obtain the excitation spectrum 
in the 7r-fiux (ground state) sector, we fix the gauge 
by choosing Ui tX = i,Ui y = (— 1) J . The corresponding 
Hamiltonian is given by 

Hq = J2i txfjfi+x + ty{—l) l fi fi+y ~ J$fifi 



which describes a p-wavc superconductor of spinlcss 
fermions [13]. Here t\ = J\ + J' x and A a = J\ — J' x . 
(Note that in Rcf. [l4j, the pairing terms are absent 
due to "projective symmetries".) In terms of the Bloch 
states, /k = ^2 i exp(—ik-r i )f i /^/N, Eq. (TTU)) in momen- 
tum space is given by 

Ho=J2 ®l H ^> $ k = (ft, fl +Q , /-k, /-k-q), (11) 



where Q = (tt, it) and the summation is over only half 
of the Brillouin zone since k is equivalent to k + Q. In 
general, the analytical form of the eigenvalues of the 4x4 
matrix Hk are complicated. Due to time reversal sym- 
metry, however, it is straightforward to derive the quasi- 
particlc excitation spectrum as follows: 
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where g±^ = (J 2 ± J' x ) cos 2 k x + (J 2 ± J' y ) sin 2 k y and 
,9k = (Jx + J' x ? cos 2 k x + (J y + J' y ) 2 sin 2 ky. 

When J 5 = 0, E +M = 4{J 2 cos 2 k x + J 2 sin 2 ky) 1 / 2 and 

E-^ = 4:(J X 2 cos 2 k x + J' y 2 sin 2 ky) 1 / 2 are the energies of 
the d and d! Majorana fermions respectively since they 
decouple from each other. Both spectra are gapless at 
nodes ±K = ±(7r/2,0), around which the spectrum is 
linear in momentum; the excitations are masslcss Dirac 
fermions. However, the massless fermions are unstable in 
the sense that additional small, local (further neighbor 
hopping and pairing) terms can gap the nodes [23| . 

When < J 5 C A„, (so the ground state lies in the 
7r-flux sector,) it is clear that E + ^ is always gapped. The 
conditions for i?_,k to have gapless excitations are: 

J x J' x cos 2 k x + JyJ' y sin 2 k y = J 5 2 /4, (13) 
(JxJ'y — JyJ' x ) cos k x sin k y = 0. (14) 

For simplicity, we consider the case in which J x , J y 3> 
J5 > so that J 5 <C A„ is satisfied for arbitrary J' x and 
J'y . From the two conditions we obtain the phase diagram 
shown in Fig. [T] as a function of J' x and J' y : (i) When 
J' x > J?/(4J X ), J'y > J?/(4J„), and J' x /J' y J x /J y , the 
fermion spectrum has eight Dirac nodes at (±7r/2±# x , 0) 
and (±n/2,±9 y ) with 9 X = arcsin( Js/^/A J\ J' x ) ; (ii) 
When J' x > J 2 /(4J X ) and J' y < J 2 /(AJ X ), there arc four 
Dirac nodes at (±tt/2 ± 9 X , 0); (iii) When J' x < Ji/{AJ X ) 
and J'y > J 2 /(AJ X ), there are also four Dirac nodes 
but at (±7r/2,±0„); (iv) When J' x < J|/(4J X ) and 
J'y < J 2 /(AJ y ), all fermionic excitations are gapped. 

For < J5 -C A„, the Dirac nodes, if they exist, are 
topologically stable in the sense of Wigner-Von Neumann 
theorem. Any weak translational and time reversal in- 
variant perturbation only shifts the positions of the nodes 
[23I ]. Consequently, the present phase with Dirac nodes 
is characteristic of a stable quan tum phase of matter, i.e. 
an algebraic spin liquid [l8, lj|, 2p| . 
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FIG. 1: The quantum phase diagram of the Gamma matrix 
model as a function of J' x and J' y in the case J5 <C J x ,J y , 
where the ground state lies in the uniform 7r-fiux sector. The 
Dirac nodes in the phase diagram are topologically stable. At 
the critical (red) line J' x /J' y = J x /J y , fermions form a Fermi 
surface (FS). The other two critical (blue) lines are defined as 
J' x = Jf /(4J») or J' y = J|/(4J„). 



in small J5 limit can be adiabatically connected to the 
large J5 gapped phase. 

Discussion: It is straightforward to generalize the 
GMM to other lattices in 2D and also to higher dimen- 
sions. For the 2D triangular lattice, a quantum spin-7/2 
(or spin-i-i-i [ljj]) model can be defined via the seven 
Gamma matrices. Due to the non-bipartiteness of the 
lattice, the model spontaneously breaks TRS. Thus, it is 
expected that a chiral spin liquid [22f could be realized 
in this model on a triangular lattice [23f . 
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It is worth noting that along the critical line J' x / J' y = 
Jx/Jy and J' x > J^/(4J X ), corresponding to the red line 
in Fig. [TJ the discrete nodes broaden into a line of nodes. 
In short, in this special case the present model realizes 
a spinon Fermi surface. Since spin-3/2 operators can be 
written as a bosonic bilinear in term of two Schwingcr 
bosons, S a = b\(jg S ,b s < /2, s =T,| with the constraint 
frj&l + frjbj = 3, the present spin model can be written as 
a bosonic model, albeit one with four-body interactions. 
It follows as a corollary that by tuning a single coupling 
constant to a critical value, a purely bosonic model can 
exhibit an emergent Fermi surface. 

Upon approach to the critical line J' x = Jf/(4J X ) 
from the eight node phase, both of the two Dirac cones 
at ±(7r/2 + 8 X ,Q) approach (ir, 0) leading, at critical- 
ity, to a single node with the unusual dispersion w 

J~Aq%. + Bqf) for small |q| = |k - (7r,0)|, where A and 

B are constants depending on J' x , J\, and J 5 . Another 
two nodes at ±(ir/2 — X ,0) approach (0,0) leading to 
the same unusual dispersion. Similar physics is obtained 
at the other critical line J' y = Jf/(4J y ). 

Large- J5 limit: The ground-states for I/J5 ~ form 
the low energy manifold Sf = ±3/2 for each i. Defin- 
ing an effective spin-1/2 S, such that Sf = ±3/2 cor- 
responds to uf = ±1, and employing degenerate per- 
turbation theory in J\ and J' x , the effective Hamilto- 
nian can be shown be H cS = JcffJ2i cr f' J i +x 'i+x+y' T t+y^ 
Jeff = (Jx - J'xYiJy - J'y) 2 7(16^|), which is exactly 
Wen's plaquette model 0]. In terms of fluxes, H c $ = 
J2 i J e s exp[i(j>i\. Consequently, the ground state is still 
in the 7r-flux sector. We expect that the gapped phase 
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